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ABSTRACT 

The Smithsonian Hectospec Lensing Survey (SHELS) is a window on the star formation 
history over the last 4 Gyr. SHELS is a spectroscopically complete survey for i? t ot < 20.3 over 
4D°. We use the 10k spectra to select a sample of pure star forming galaxies based on their 
Ha emission line. We use the spectroscopy to determine extinction corrections for individual 
galaxies and to remove active galaxies in order to reduce systematic uncertainties. We use the 
large volume of SHELS with the depth of a narrowband survey for Ha galaxies at z ~ 0.24 
to make a combined determination of the Ha luminosity function at z ~ 0.24. The large area 
covered by SHELS yields a survey volume big enough to determine the bright end of the Ha 
luminosity function from redshift 0.100 to 0.377 for an assumed fixed faint-end slope a = —1.20. 
The bright end evolves: the characteristic luminosity L* increases by 0.84 dex over this redshift 
range. Similarly, the star formation density increases by 0.11 dex. The fraction of galaxies with 
a close neighbor increases by a factor of 2 — 5 for L^ a > L* in each of the redshift bins. We 
conclude that triggered star formation is an important influence for star forming galaxies with 
Ha emission. 
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Determining the star formation history of the 
Universe is a crucial part of understanding the for- 
mation and evolution of galaxies. Exploration of 
the global star formation history has two compo- 
nents: (i) measurement of the star formation den- 
sity over time and (ii) understanding the physical 
processes that drive star formation. Here we use 
a large, moderate-depth spectroscopic survey to 
address both issues: (?) wc determine the star for- 
mation density over the last 4 Gyr using the Ha 
emission line as star formation indicator and (ii) 
we investigate the possible influence of galaxy in- 
teractions on the Ha luminosity function. 

There is abundant observational evidence for 



an order of magnitude increase in the star forma 
tion d ensity since redshif t z ~ 1 — 2 (jLillv et al 



1996;_ Madau et all Il996t and the compi lations 



of iHopkins! 12004 iHopkins fc Beacomll2006l ). Ma- 
jor mergers, tidal interactions, gas removal from 
conversion into stars, and/or ram pressure strip- 
ping may explain the decrease in the star forma- 
tion. The challenge is deciding which of these pro- 
cesse s are important in quenching of star forma- 



tion (jBell et alj|2005h 



The decline in star formation density coincides 
with a rapid decrease in the characteristic lumi- 
nosit y of galaxies (L* ) in the rest-frame £/-band 



(e.g. Illbert et all 120051 : IPrescott et all l2009h . A 



1 Smithsonian Astrophysical Observatory, 60 Garden 
Street, Cambridge, MA 02138, USA 

2 Brown University, Department of Physics, Box 1843, 
Providence, RI 02912, USA 



decrease in the number of merging systems can 
explain the decrease of the c h aracteristic luminos - 
ity L* (|Le Fevre et alJl200dh . ISobral et all (|2009h 
find a strong morphology-Ha luminosity relation 
for mergers and non-mergers. The characteristic 
luminosity L* defines a critical switch-over lumi- 
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nosity between the mergers and non-mergers; the 
mergers are more luminous. 

Studies of close pairs show that enhancement 
in the star formation rate is large st for galax- 
ies in major pairs (|Am| < 0.7 — 2 



Wood s et al.l 
Ellison et al.l 1200 



200a IWoods & Gellerl 12007 
and that the average star formation rate in a 
galaxy increases with decreasing projected sepa- 
ration (|Li et al.ll2008l ). Simulations of interacting 
and merging galaxies reveal that the interactions 
can trigger short powerful bursts of star formation 
by forcing subst antial fractions of the gas into the 
central regions ( Mihos fc Hernquist 19961 ) . 

Systematic effects dominate the comparison of 
star formation rates determined from different 
star formation indicators like the rest-frame ultra- 
violet (UV) and Ha. Hence, to study the varia- 
tion of the star formation density with time, the 
use of a single star formation indicator is best. 
The rest-frame UV spectrum of a galaxy directly 
measures the population of newbor n stars (e.g. 



Lilly et all Il996l : iTrever et all Il998f ). However 



the rest-frame UV is strongly attenuat ed (e.g. 



Cardelli et all 119891 ICalzetti et all l2000h . The 



most-direct optical indicator is the Ha emission 
line emitted by gas surrounding the e mbedded star 
forming region (e.g. lKennicuttlll998l ). The Ha line 
is also affected by attenuation-albeit less than the 
UV-which can be corrected using spectroscopy. 
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2009) to determine the Ha luminosity function 
parameters over a range of redshifts. Despite 
the depth of the narrowband surveys, measure- 
ments of individual luminosity function parame- 
ters and the star formation density are not well- 
constrained. 

Narrowband surveys lack spectroscopy for the 
faint Ha emitting galaxies. Thus, general as- 
sumptions about stellar absorption, extinction 
corrections, contributions by active galactic nu- 
clei (AGNs), or interloper contamination need to 
be made for the sample as a whole rather than 
for each galaxy. These issues may lead tosys- 



tematic uncertainties. iMassarotti et alJ ((2001) 



tion introduces a systematic underestimate of the 
extinction-corrected star formation density. A 
spectroscopic survey does not suffer these limita- 
tions, although it is usually limited in its depth. 

Several spectroscopic Ha surveys exist (e.g. 

Gallego et al.ll995t|Tresse fc Maddoxll998l:ISullivan et al 



2000 1 : iTresse et alJ 120021: IPerez-Gonzalez et al. 
20031 IShim et al.ll2009h . Both iGallego et alJ and 



IPerez-Gonzalez et al.l use the Universidad Com- 
plutense de Madrid (UCM) survey. This sur- 
vey covers an extremely wide area on the sky 
(472 D ). However, it is limited to a very low 
redshift (z ma ,v ~ 0. 045). For their Ha survey, 
ISullivan et all (|2000h use galaxies selected from 
UV imaging in a 2.2 D° field. The other surveys 
have an area < 0.25 D°. Thus, most surveys are 
too limited in volume to overcome cosmic vari- 
ance. 

The Smithsonian Hectospec Lensing Survey 
(SHELS) is a spectroscopic survey covering 4D° 
on the sky to a limiting i? -band magnitude 



R tot = 20.3 (|Geller et al.ll2005h . We use SHELS 



show that applying an average extinction correc- 



to obtain a consistent determination of the star 
formation history over the last 4 Gyr based on the 
Ha emission line over a relatively large area and 
redshift range. 

The spectroscopy enables us to reduce system- 
atic uncertainties by allowing an individual galaxy 
extinction correction. We can also remove individ- 
ual AGNs rather than applying a global correc- 
tion factor for contamination by AGNs as is done 
in narrowband surveys. We use the large survey 
area to determine the characteristic luminosity L* 
of the Ha luminosity function and associated sys- 
tematic uncertainties. 

We discuss the SHELS spectroscopic data in 
Section [2] In Section [3] we introduce our Ha 
sample selection. We combine our i?-band se- 
lected Ha sample with th e narrowband Ha sur- 
vey of IShiova et al. ( 20081 ) in Section 0] to obtain 
a jointly-determined Ha luminosity function at 
z ~ 0.24. Sections [5] and [6] discuss the deriva- 
tion and evolution of the luminosity function and 
star formation density, respectively, over the past 
4 Gyr. We include an investigation of the influence 
of our selection criteria on the derivation of the 
luminosity function parameters. In Section [7] we 
examine the stellar age of the star forming galax- 
ies and the influence of galaxy-galaxy interactions 
on these galaxies. We summarize our results in 
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Section [3 

Throughout this paper we assume a flat Uni- 
verse with H = 71 kms -1 Mpc -1 , Q m = 0.27 and 
f^A = 0.73. All quoted magnitudes are on the AB- 
system and luminosities are in ergs -1 . 

2. SHELS observations 

We constructed the SHELS galaxy catalog from 
the i?-band so urce list for the F2 fie l d of th e Deep 
Lens Survey (jWittman et alj 120021 120061) . The 
DLS is an NOAO key program covering 20 D° in 
five separate fields; the 4.2 D° F2 field is centered 
at a = 09' l 19 m 32.4 s and S = +30°00'00". We ex- 
clude regions around bright stars (~ 5 % of the to- 
tal survey) resulting in an effective area of 4.0 0°. 
We use surface brightness and magnitude to sep- 
arate stars from galaxies. This selection removes 
some AGN. 

Photometric observa tions of F2 were ma de with 
the MOSAIC I imager (jMuller et al.lll998f) on the 
KPNO Mayall 4 m telescope between 1999 Novem- 
ber and 2004 November. The i?-band exposures, 
all taken in seeing < 0.9" FWHM, are the ba- 
sis for the SHELS survey. The effective exposure 
time is about 14,500 seconds and the 1 a surface 
brightness limit in R is 28.7 magn itudes per square 
arcsecond. IWittman et al.l (|2006l ) describe the re- 
duction pipeline. 

We acquired spectra for the g alaxies with the 



emiss ion-line objects (see also iFabricant et al 



2005) 



Hectospec f iber-fed spectrograph (jFabricant et al 
19981 120051 ) on the MMT from 2004 April 13 
to 2007 April 20. The spectrograph is fed by 
300 fibers that can be positioned over a 1° field. 
Roughly 30 fibers per exposure are used to deter- 
mine the sky. The Hectosp ec observation planning 
software ( Roll et al.lll998l ) enables efficient acqui- 



sition of a magnitude limited sample. 

The SHELS spectra cover the wavelength range 
A = 3,500 - 10, 000 A with a resolution of -6 A. 
Exposure times ranged from 0.75 hours to 2 hours 
for the lowest surface brightness objects in the 
survey. We reduced th e data with the stan- 
dard Hectospec pipeline dMink et al. 2007) and 
derived redshifts with RVSAO (Kurtz & Mink 



1998 ) with temp l ates c onstructed for this purpose 
([Fabricant et alj 120051 ). We have 1,468 objects 
that have been observed twice. These repeat ob- 
servations imply a mean internal error of 56 km s _1 
for absorption-line objects and 21kms~ 1 for 



Fabricant et alJ (|2008l ) describe the technique 



we use for photometric calibration of the Hec- 
tospec spectra based on the particularly stable in- 
strument response. For galaxies in common be- 
tween SHELS and SDSS, the normalized Ha line 
fluxes agree well in spite of the difference in fiber 
diameters for the Hectospec (l'/5) and the SDSS 
(3"). For high-signal-to-noise SHELS spectra, the 
typical uncertainties in emission line fluxes are 
18%. 

SHELS includes 9,825 galaxies to the limiting 
apparent magnitude. The overall completeness of 
the redshift survey to a totafj] i?-band magnitude 
of i?tot < 20.3 is 97.7%, i.e. 9,595 galaxies have a 
redshift measured; the differential completeness at 
the limiting magnitude is 94.6%. The 230 objects 
without redshifts are low surface and/or faint ob- 
jects, or objects near the survey corners and edges. 
M. J. Kurtz et al. (2010; in preparation) includes 
a detailed description of the full redshift survey. 

The SHELS survey also includes 1,852 galaxies 
with 20.3 < R < 20.6, for which we have measured 
a redshift; the total sample of galaxies with 20.3 < 
R < 20.6 is 3,590, i.e. the survey is 52 % complete 
in this magnitude interval. The completeness is 
patchy across the field. 

The F2 field contains an atypical under-dense 
region at the lowest redshifts because the DLS 
fields are selected against nearby clusters at z < 
0.1. We show the redshift distribution of our Ha 
galaxies in Figure [IJ 

2.1. The i?-band k + e-corrections 

To calculate the absolute i?-band magni- 
tude Mr we determine the appropriate k + e- 
corrections. The k + e-correction converts the 
observed absolute magnitude to the rest-frame of 
the galaxy, correcting for redshift and evolution. 
We use the k + e-corrections for 9 types of galax- 
ies: bright cluster (BCG), elliptical (E), SO, Sa, 
Sb, Sbc, Sc, Sd and irregular (Irr) galaxies deter- 
mined by J. Annis (priv. comm.jj. We use the 



1 The total magnitude is the SExtractor dBertin fc Arnoutsl 
1996) MAG_AUTO as opposed to an aperture magnitude. 

2 The table with the corrections for the SDSS filter set as 
function of galaxy type and redshift can be obtained from 
http: //home . fnal .gov/- annis/astrophys/kcorr/kcorr .html 
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Fig. 1. — Redshift cone diagram for the galaxies in the final sample: i?tot < 20.3, S/Nhc* > 5 and /hq > 
2Q-15.5 er g S -i cm -2 AGNs have not been removed from this sample. The large-scale structure is apparent 
with extended low-density regions and well-populated narrow structures. 



corrections for the SDSS r'-filter as a function of 
redshift and (g' — i')-color because the SDSS r'- 
filter is similar to the i?-filter used for the DLS. 
We obtain {g' — i') by cross- matching our catalo g 
with SDSS DR6 (|Adelman-McCarthv et aT1l2008fh 
For those galaxies not found in SDSS DR6 (these 
galaxies are either unresolved or below the surface 
brightness limit in SDSS) we convert the (V — R) 
from 61 galaxies in the DLS to (g' — i'). For 42 
galaxies we cannot determine or derive {g 1 — i') 
due to the proximity of another object; we assume 
that these are Sa galaxies. We interpolate the 
models in redshift to obtain the k + e-corrections 
for each galaxy type determined by its (<?' — i') 
and redshift. 

3. Ha sample selection 

We use SHELS to construct Ha luminosity 
functions over the redshift range 0.010 < z < 
0.377 (Table Q]). Here, we describe the determina- 
tion of our final emission-line luminosities and the 



discrimination between pure star-forming galaxies 
and AGNs. 

3.1. Emission-line measurements 

The emission-line flux emanating from star- 
forming regions is affected by the absorption-line 
spectrum from the underlying stellar population. 
The absorption mostly affects the measurements 
of the hydrogen Balmer lines. To measure the 
emission-line flux we thus remove the contribution 
of the stellar population. 

(|2004[ l contin- 



We use the Trcmonti et al 



uum subtraction method to correct for the stel- 
lar absorption rather than applying a con stant, 



global correctio n (e.g., Hopkins et al. 20031 ). The 



iTremonti et alJ method removes the stellar contin- 
uum by fitting a linear combination of template 
spectra rcsampled to the correct velocity disper- 
sion. The method also accounts for redshift and 
reddening. The template spectra are based on 
single stellar population models generated by the 
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Fig. 2. — Fraction of light contained in the l'/5 
fibers as a function of redshift. We indicate the 
fraction of galaxies with more than 20% of the 
light contained in the fiber (horizontal dashed 
line) for each redshift bin (vertical dashed lines) 
we use to construct the Ha luminosity functions. 
The galaxies have i? to t < 20.3, S/N Hq > 5 and 
/h„ > 10~ 15 - 5 ergs -1 cm -2 . The AGNs have been 
removed. 
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Fig. 3. — Attenuation at Ha as function of 
observed Ha luminosity. The black line indi- 
cates the least-absolute-deviates fit to the gray 
points. We indicate A^ a — (solid hor- 
izontal line) and a commonly assumed value 



of A 



e.g. iTresse fc Maddoxl Il998t iFuiita et all 120031 : 
Lv et alJl2007tlSobral et alj|2009l) 



1, dashed horizontal line: 



popul ation synthesis code of Bruzual fc Charlotl 
(|2003l ). We use models with 10 different ages 
(0.005, 0.025, 0.1, 0.3, 0.6, 0.9, 1.4, 2.5, 5 and 10 
Gyr) at solar metallicity. 

We determine the emission-line fluxes from the 
continuum-subtracted spectra by integrating the 
line flux within a top-hat filter centered on the 
emission-line. We remove any local over- or under- 
subtraction of the continuum by subtracting the 
mean of the flux-density at both sides of the filter. 
Next, we determine the continuum level by tak- 
ing the mean of the flux-density at wavelengths 
bluer and redder than the emission-line on the 
best-fit continuum model. Finally, we determine 
the absorption contribution of the underlying stel- 
lar population using the same top-hat filter but on 
the best-fit model; we remove the flux contributed 
by the continuum. 

The Hectospec fibers have a fixed diameter of 
l'/5. At all redshifts where Ha is observable (and 
in particular at the lowest redshifts) the fiber does 
not cover the entire galaxy. Hence, we use an aper- 
ture correction 

j\ — ^Q-0.4(m tataJ -m aber ) 

to correct for the fiber-covering fraction. Figure [2] 



shows the fraction of light, I /A, contained in the 
fiber as a function of redshift. 



Kewlev et al 



(120051 ) show that a spectrum 
measuring at least 20 % of the galaxy light avoids 
substantial scatter between the nuclear and inte- 
grated SFR measurements. The overall major- 
ity of galaxies from SHELS have a light-fraction 
1 1 A > 20% (Figur e [2]). 

iFabricant et al.l (2008) compared the Ha and 
[On] emission-line fluxes from SHELS with SDSS 
DR6 after making an aperture correction. They 
found excellent agreement between the two sur- 
veys, even though the fibers of the SDSS spec- 
trograph are 3" in diameter. Moreover, most of 
the SDSS galaxies are at low redshift (z < 0.14) 
where we have the largest fraction of galaxies with 
a light-fraction less than 20%. 

There is no dependence of final Ha luminosity 
on the light-fraction. We are thus confident that 
the use of these aperture corrections does not af- 
fect the final results even when the covering frac- 
tion is small. 
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Fig. 4. — Ha luminosity (corrected for extinc- 
tion) as a function of redshift. The solid line in- 
dicates the additional selection criterion /hq > 
10 -155 ergs -1 cm' 2 . Vertical dashed lines show 
the edges of the redshift bins used to construct 
the SHELS Ha luminosity functions. 



3.2. Extinction correction 

The light from star forming regions in a galaxy 
is often heavily attenuated. To determine the in- 
trinsic SFR of a galaxy we must remove the effects 
of attenuation. 

We calculate the attenuation by comparing the 
observed value of the Balmer decrement (corrected 
for stellar absorption) with the theoretical value 
(iW/H/3 = 2.87 for T = 10 4 K and case B re- 
combination; Table 2 of ICalzetti|[200lh . The in- 
trinsic flux is /intr(A) = /ob s (A)10 a4A \ where A x 
is the wavelength-dependent extinction. A\ is 



Ay 



= *(A) 



k(X)E(B - V) gas 
2.51ogi? a/ 3 



fe(H/3) - fc(Ha) 



(2) 



where R a p is the ratio of the attenuated-to- 
intrinsic Balmer line ratios, k(H/3) — fc(Ha) is 
the differential extinction between the wave- 
lengths of H/3 and Ha, and fc(A) is the extinc- 
tion a t wavelength A. We apply the lCalzetti et al 



(|2000l ) extinction law, which has k(V) = 4.05, 
fc(Ha) = 3.325 and fc(H/3) = 4.596. 

Figure [3] shows the attenuation as a function of 
observed Ha luminosity for galaxies with both Ha 
and H/3 at a S/N > 5. We use the relation between 
An a and Lhq as determined from a least-absolute- 
deviates fit to the high-S/N data-points with ob- 



served luminosities 40.5 < log I/Ha < 41.5 (gray 
points) for galaxies where S/N H/3 < 5 or where 
the observed equivalent width of H/3, OEWh/3 < 
1 A (uncorrected for stellar absorption) . We limit 
OEWh/3 to avoid galaxies with excessively large 
attenuation resulting from a very small (noise- 
dominated) H/3 flux compared to Ha. We assume 
that galaxies with An a < to have no attenuation 
and assign Ab_ u = to these galaxies. 

3.3. Sample definition 

Figure U shows the Ha luminosity as a function 
of redshift. Below /n a = 10 -15 5 ergs -1 cm -2 the 
number of galaxies decreases rapidly. We impose a 
constant Ha flux limit /h q > 10 -15,5 ergs -1 cm -2 
on the sample (after corrections for stellar absorp- 
tion and attenuation) because we are only com- 
plete to this flux. We apply this criterion in addi- 
tion to the magnitude limit (i?tot < 20.3) and the 
S/N Ha > 5 requirement. 

3.4. AGN classification 

The presence of an active nucleus in a galaxy 
contributes to the (app arent) star-formation in the 
galaxy. For example, iPascual et al. (|200lh find 
that approximately 1 5 % of the luminosity density 
of the UCM survey (|Gallego et al.1119951) results 



from g alaxies identified as AGN. Westra fc Jonesl 
(2008) find a 5 % contribution for their survey. 

We use the demarca t ions o f pure star formation 
from Kauffm ann et al. ( 20031) a nd of extreme star- 
burst from lKewlev et alj (|2001[ ) to identify galax- 
ies as pure star forming, AGN, or a combination 
(composite galaxies) based on the line ratios of 
[O in] /H/3 and [Nil] /Ha (Figure [5}. 

We select all galaxies with [O in] A5007 and 
[N n] A6585 detected with a signal-to-noise ratio 
(S/N) > 3. If the galaxies have both Ha and H/3 
detected with S/N > 3, we identify them as pure 
star formi ng galaxies when their line ratios are 
below the iKauffmann et alj re lation, pure A GNs 



rela- 



when the ratios are above the iKewlev et al. 
tion, and composites when they lie between the 
relations. 

For galaxies with either Ha or H/3 undetected 
(S/N < 3), we use the 3 a value for the line flux 
to calculate the line ratios. These ratios are lower 
limits. We classify these galaxies as composite or 
AGN; some of the composite galaxies might be 
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Fig. 5.— BPT (after lBaldwin et alJl98lf > diagram 
for SHELS. The solid blue and dashed red lines in- 
dicate the demarcation o f pure star formation from 




Kauffmann et al.l ( 20031) and of extreme starbursts 



( 20011 ) . We classify the galaxies 



as: pure star forming galaxies (black diamonds), 
AGNs (red squares), and composite galaxies (blue 
triangles). We indicate galaxies or AGNs with ei- 
ther Ha or H/3 undetected, i.e. S/N < 3, as lower 
limits (green arrows). 



AGNs. 

We identify a separate class of broad-line AGN. 
The width of these broad Balmer- lines extends be- 
yond the limited-width top-hat filter used for mea- 
suring the line fluxes (Section 13. In some cases 
the [N n] AA6550,6585 lines are not distinguishable 
from the Ha line in a spectrum with a very broad 
Ha line. 

Inspecting each spectrum would be time- 
consuming. Hence, we fit the Ha and H/3 lines 
in the continuum-subtracted spectra in an au- 
tomated way and individually inspected each 
candidate broad-line AGN. We fit both lines si- 
multaneously with the assumption that the full- 
width-half-maximum (FWHM) of the line profile 
is the same for both lines. Candidate broad- 
line AGNs have a peak of both Ha and H/3 
> 5 x 10 -18 ergs -1 cm~ 2 A -1 above the contin- 
uum residuals (which avoids the inclusion of noise 
peaks) and a FWHM of the Gaussian component 
of the line profile (we use a Gaussian convolved 
with the instrumental profile as our line profile) 
before convolution larger than 14 A. From these 
candidates, we select the galaxies that are genuine 



broad- line AGNs. 

The fraction of galaxies identified as AGN 
and/or composite over the redshift ranges 0.010- 
0.100, 0.100-0.200, 0.200-0.300 and 0.300-0.377 for 
an Ha luminosity limited sample (logL > 41.18; 
lowest Ha luminosity at z = 0.377) is 5.9, 6.6, 5.3 
and 5.2%, respectively. 

The fraction of AGN is more or less constant 
with redshift. However, we cannot draw any con- 
clusions about the evolution of the AGN-fraction 
as a function of redshift. We removed stellar ob- 
jects from the initial sample and thus may have 
inadvertently removed AGNs particular at greater 
redshifts. 



4. 



The Ha luminosity function at redshift 
- 0.24 



on 



The recent advent of wide-field cameras 
telescopes has aided searches for star forming 
galaxies by increasing the area (an d hence vol- 
ume) of narrowb a nd sur veys, e.g. iFuiita et al 



Westra k, Joned ( 20081 ). and many 



(120031). iLv et al.( (120071) . IShiova et all (I2008D . 

more. This 



technique ha s recently been exte nded to the near- 
infrared, e.g. ISobral et al. I (|2009h . 

A narrowband survey efficiently probes the 
faint end of the luminosity function which is hard 
to explore in a spectroscopic survey. In contrast, a 
spectroscopic survey can cover a larger volume and 
sample the rare luminous galaxies at the bright 
end of the luminosity function. 

Here, we combine the strength of a narrowband 
survey-the ability to go deep-with that of our 
broadband selected spectroscopic survey-coverage 
of a large volume-to determine a well-constrained 
luminosity function at z ~ 0.24. For the narrow- 
band suryev__weuse_the publicly avail able data 
from IShiova et all (|2008l hereafter Isoj ) together 
with t hat of the Cosmic E volution Survey (COS- 
MOS^; ICapak et aT]|2007t ) which formed the basis 



of the survey of lS08l . We use the spectroscopic sur- 
vey of SHELS for the bright end of the luminosity 
function. 



The IS08I and SHELS surveys use different ap- 
proaches (imaging versus spectroscopy). A com- 
parison of the data allows a consistency check of 



3 The COSMOS catalog can be downloaded from 
http : //irsa . ipac . caltech . edu/data/C0SM0S/tables/cosmos_phot_20060103 . tb] 
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Fig. 6. — Comparison of SHELS with ISOa for galaxies with the OEW HQ+ p^ n j > 12 A not corrected for 

[N n] for both surveys. The data are: IS08l (solid black circles), SHELS in the redshift range of S08 (open red 
squares), and SHELS at 0.01 < z < 0.15 (open blue triangles). The Ha luminosity, Lhq, for both surveys is 
corrected for [N II] . 



the aperture corrections applied to the SHELS 
data. 

We construct the Ha lu minosity function over 
the redshift range of [sol (0.233 < z < 0.251) 
based on the catalog with emission-line fluxes de- 
termined in Section 13.11 (which already include 
corrections for underlying stellar absorption), red- 
shifts, extinction corrections from Section T3. 21 and 
removal of composites and AGNs (Section 13.41) . 
Constraining SHELS to the same redshift range 
yields a sample of 192 SHELS galaxies at 0.233 < 
z < 0.251. 

4.1. Data comparison 

Figure [6] shows the Ha luminosity, Ha rest- 
frame equivalent width (REW), and t he 3" aper- 
ture absolute i?-band magnitude from |S08| (solid 
black circles) matched to the selection criteria of 
SHELS, Rtot < 20.3. We also show the SHELS 
data (o ven red squares) for the redshift ra nge c ov- 
ered bvlsoil (0.233 < z < 0.251). To match|solwe 



require an observed equivalent width of Ha com- 
bined with [N II] > 12 A (as per the selection crite- 
ria oflSQSlL and f Ha > 10~ 15 - 5 ergs" 1 cm- 2 . The 
strengths of both surveys are immediately appar- 
ent. SHELS includes the highest luminous galax- 
ies, IS08l probes the faint end of the luminosity 
function. 

S08l measure magnitudes from a 3" aperture, 



scale them to the total i'-band magnitude, and cal- 
culate Ha fluxes. We measure fluxes in a similar 
way. We use spectra taken with a l'.'5-fiber aper- 
ture scaled to the total .R-band magnitude. The 
data from the lower redshift range 0.01 < z < 0.15 
of SHELS (Figure[5J open blue triangles) show that 
the relation between the Ha luminosity and the 3" 
total .R-band magnitude of the two surveys is sim- 
ilar. The scaling of the Ha flux from the limited- 
aperture magnitude to the total magnitude intro- 
duces no systematic biases and is consistent with 

Isoi . 

If we constrain the data of |S08| to fn a = 



10 



-15.5 



erg s cm , a difference at the bright end 
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T 42.5 




-21 

M (3" aperture) 



Fig. 7. — Comparison of the total Ha luminosity 
with the 3" aperture i?-band magnitud e of S HELS 
(open red squares) and the survey of IS08l (solid 
black circles). The galaxies have OEW Hq+ |-j^ n j > 

12A, a flux-limit of f Ka > 1CT 15 - 5 ergs -1 cm -2 , 
Rtot < 20.3, and 0.233 < z < 0.251. We show 
the median Ha luminosity for 0.5 magnitude- wide 
bins for SHELS (red line) and |S08| (black line). 
The SHELS galaxies shift toward greater Lh q at 
fixed Mr 

becomes apparent (see Figure [7]). T here are more- 
luminous galaxies in SHELS than in S08. This dif- 
ference results from two effects: (i) SHELS probes 
a larger volume, and (ii) the Ha fluxes determined 
from narrowband surveys can easily underestimate 
the true line flux. Galaxies with redshifts that 
place the Ha line in the wings of the filter under- 
estimate the mean recovered Ha flux. 

To exa mine the redshift distribution of the 
galaxies in |S08| . Figure [5] shows the red shift distri- 
butio n of the 10k zCOSMOS catalog (|Lillv et al 



20071) in combination with the filter transmission 



curve of the NB816 normalized to the maximum 
throughput. Any galaxy with a redshift placing 
it in the wings of the narrowband filter has its 
Ha flux underestimated far more than the 21 % 
|S08| use to correct their line fluxes. In the COS- 
MOS field the galaxies tend to be at redshifts to- 
wards the red edge of the filter. In this case, the 
[N 11] A6585 line (the strongest of the two [N 11] 




40 3 



Fig. 8. — Redshi ft distribution of the galaxies in 
zCOSMOS DR2 (|Lillv et al.ll2007l . histogram), the 
rela tive NB816 filter transmission curve used by 
|S08| (gray solid line), and where the transmission 
of the NB816 filter is 50% of its maximum (dotted 
lines). 



lines that straddle Ha) barely contributes to the 
flux probed by the filter. Both the underestima- 
tion of the Ha flux and over-correction for [N 11] 
can explain the difference in the distribution of Ha 
fluxes in Figure [7J 

Despite this differ ence , we can still use the 
fainter galaxies from |S08| to determine the faint- 
end slope of the Ha luminosity function. The sys- 
tematic underestimation of fluxes causes a shift in 
the luminosity function which affects the determi- 
nation of the characteristic luminosity (i.e. bright 
end), not the faint-end slope. 

4.2. Derivation and fit 



We fit a Schechter function ( Schechtei 19761 ) to 
the SHELS andlS08| data. The Schechter function 



is 



(j)(L)dL = c/)* 



exp 



L 

L* 



L 

L* 

(3) 

where a is slope of the faint-end part, L* is a char- 
acteristic luminosity, and <jf is the normalization. 
Throughout this paper the units for the Schechter 



parameters L* and <f>* are ergs and Mpc , re- 
spectively, a is dimensionless. 

To combine the two data sets we use th e non- 
parametric 1/V max method (|Schmidtlll968l ) to de- 
http://www.naoj.org/Qbserving/instruments/scam/txt/NB8i6.txttermine the Schechter parameters. The number 



4 zCOSMOS DR2, which can be obtained from the ESO 
archives. 

5 The filter profile is available at 
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Log I_ h (ergs s ') -> 

Fig. 9. — The 1/Vmax data-points for 
SHELS (solid squares), SHELS without the 
OEW HQ+ jj^ n j > 12 A criterion (solid dia- 
monds), and S08 (open squares). The galaxies 
have OEW Hq+ [ Nh ] > 12 A, R to t < 20.3, and 

0.233 < z < 0.251. The thick solid line indi- 
cates the combined fit of SHELS and S08 with 

Cmb = -3- 05 ± °- 09 > a = ~ 1A1 ± °- 03 ' and 
logL* = 42.14 ± 0.08. The thin lines indicate 

the luminosity function for the SHELS and |S08| 

data separately with </>g HELS = — 3.11 ± 0.09 

(dashed line) and <^g 08 = —2.91 ± 0.09 (dotted 

line), respectively. Both luminosity functions also 

have a = -1.41 ±0.03 and logL* = 42.14 ±0.08. 

density of galaxies for each luminosity bin j with 
a width of A log L is 



4>(L j )b\DgL=Y, 



(4) 



where W(x) = 1 when the luminosity is enclosed 
by bin j and W(x) — otherwise, and Vi is the 
volume sampled by galaxy i. The uncertainties in 
the bins are Poisson errors 



E 



W(U) 
Vi 



(5) 



Figure [9] shows the luminosity function for the 
combined data set (thick solid line). For both sur- 
veys we apply the selection criteria i?tot < 20.3, 



OEW 



> 12 A, and 0.233 < z < 0.251. 



Ha+[N II] 

We determine the data-points using 1/V max 
where the uncertainties are Poisson errors for both 



SHELS (solid squares) and lS08l (open squares). We 
fit a Schechter function with comm on L* and a 
to the combined SHELS and IS08l dataset. For 
this fit, we use the data with log Lnq > 41.4 
for SHELS and logi H a > 39.6 for [§M We re- 
cover a single value for a and L* of the joint 
fit: a = -1.41 ± 0.03 and logL* = 42.14 ± 0.08. 
For those values, the normalizat ion f or SHELS is 
log^SHELS = — 3.11±0.09 and for lS08l is log 6% oz = 



"SHELS 

-2.91 ± 0.09. We combine 



y S08 



and 



^SHELS 



ing a volume- weighted average. Thus, log0* omb = 
log[(1.5 x 10~ 2 - 91 +4.0 x 10" 311 )/5.5] = -3.05. 
This method of combining is one way to account 
for cosmic variance. We adopt a — —1.41 ± 0.03, 
logL* = 42.14 ±0.08, log0* omb = -3.05 ± 0.09 for 
comparison to other surveys. 

5. The Ha luminosity functions from 
SHELS 

Here we use SHELS to determine the Ha lu- 
minosity functions as a function of redshift. We 
can identify Ha in our spectra up to a redshift 
of z m ax — 0.377. We next examine the influence 
of the i?-magnitude limited survey on the deriva- 
tion of the Ha luminosity function. We limit our 
sample to /h q > 10~ 15 ' 5 ergs -1 cm -2 (see Sec- 
tion 13. 3p . Table Q] lists the number of galaxies 
satisfying each selection criterion. Figure |4] shows 
the distribution of the Ha luminosity as a function 
of redshift and the redshift bins used to construct 
the Ha luminosity functions. 

5.1. Derivation and fit 



He re we use the STY-method (jSandage et al 
1979), a parametric estimation method, to deter- 
mine the three Schechter parameters for each red- 
shift bin. The STY-method identifies the luminos- 
ity function parameters that maximize the proba- 
bility of obtaining the observed sample. The prob- 
ability V is 



N,. 



<KLi) 



(6) 



where L(z)n m ^ is the faintest luminosity where 
galaxy i at redshift zi is observable. We use a 
truncated-Newton method to maximize the natu- 
ral logarithm of V . 

Table [2] lists the fit-parameters, and Figure fTOl 
shows the results. We fit for a, L* , and 4>* (dashed 



10 



Table 1 

Number of galaxies satisfying each selection criterion. 





Criterion 


0.01 < z < 0.10 


0.10 < z < 0.20 


0.20 < z < 0.30 


0.30 < z < 0.38 


Total 




i?tot < 20.3 


461 


1949 


2746 


2114 


7270 


(it) 


S/N HQ > 5 


420 


1640 


1857 


1275 


5192 


(Hi) 


log HQ > -15.5 


369 


1441 


1702 


1186 


4698 


(iv) 


pure star forming 


322 


1127 


1268 


848 


3565 



Note. — Successive lines are a subset of the line above. We apply the selection criteria sequentially in the 
order of the Table. 



lines). We also use a fixed a = —1.20 (solid lines). 
This fixed value represents the slope over the red- 
shift range 0.05 < z < 0.20 for SHELS. This range 
has a large enough volume to sample the bright 
end of the luminosity function, while still having 
galaxies faint enough to determine the faint end 
slope. We do not consider this large redshift range 
in our further analysis. 

Narrowband surveys apply a correction to 
the total lumi nosity density for galaxies hosting 
an AGN (e.g. iFuiita et all 120031; Lv et all [2007 



Shiova et alfeoOSUWestra fc Jonesll2008h . A large 
fraction of these surveys have little or no spec- 
troscopy. Thus there is no way to separate AGNs 
from star forming galaxies. SHELS enables a di- 
rect separation (see Section |3T|) . We derive the 
Schechter parameters for the Ha emitting galax- 
ies with the AGNs removed (Table [2]). Removal of 
AGNs moves L* slightly fainter and reduces the 
normalization because AGNs are systematically 
in more luminous galaxies. Failure to account for 
this bias introduces a systematic offset. 

We determine the final uncertainties in the 
Schechter parameters by constructing 1,000 sets of 
Ha luminosities. We simulate the Ha luminosities 
by converting the observed absolute magnitudes 
into Ha luminosities using the distribution of Ln a 
as a function of Mr for each redshift bin from 
Figure 111! We redetermine the Schechter param- 
eters for each simulation using the STY-method. 
The 1 a spread in the redetermined parameters is 
the final uncertainty which includes the formal fit- 
ting uncertainty, uncertainties resulting from the 
size of the sampled volume, and the uncertainties 




□ .DID < i < 0.100 O 
0.1 DO < t < 0.200 O 
0.300 < i < 0.300 O 
0.300 < I < 0.377 O 



Fig. 11. — The distributions of logLHa as func- 
tion of Mr for each redshift bin. We use these 
distributions to assign an Ha luminosity to the 
observed absolute magnitude and to assign an ab- 
solute magnitude to a simulated Ha luminosity. 
This procedure enables us to determine the final 
uncertainties in the Schechter parameters and to 
study the influence of the survey selection criterion 
Rtot < 20.3. 



in the observed Ha luminosity. Table [5] lists the 
Schechter parameters and their uncertainties. 

5.2. Parameter evolution and impact of se- 
lection criteria 

Figure [12] compares the Schechter function pa- 
rameters of SHELS with fixed a = —1.20 with 
other Ha surveys. Evolution in the characteristic 
luminosity L* is clearly visible. 
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Table 2 

Parameters for the Ha luminosity functions. 







fixed a unconstrained a 








redshift range 


Ot 


logL* log0* a logL* 


log^ 










Pure star forming sample 








0.010 < z < 0.100 a 


— 1.20 


41.72 ±0.10 -2.86 ±0.04 -1.22 ±0.06 41.74 ±0.13 


-2.90 ± 





10 


0.100 < z < 0.200 


— 1.20 


42.09 ± 0.04 -2.97 ±0.02 -0.87 ± 0.05 41.79 ±0.06 


-2.58 ± 





05 


0.200 < z < 0.300 


-1.20 


42.52 ± 0.04 -3.29 ±0.01 -0.71 ± 0.07 42.13 ±0.06 


-2.79 ± 





05 


0.300 < z < 0.377 


-1.20 


42.83 ± 0.03 -3.59 ±0.01 -0.50 ± 0.06 42.30 ±0.05 


-2.96 ± 





04 


0.233 < z < 0.251 b 




-1.41 ±0.03 42.14 ±0.08 


-3.05 ± 





09 




Pure 


star forming sample including composites and AGNs 








0.010 < z < 0.100 a 


-1.20 


41.76 ±0.10 -2.82 ±0.04 -1.25 ±0.06 41.85 ±0.13 


-2.93 ± 





10 


0.100 < z < 0.200 


-1.20 


42.13 ± 0.04 -2.88 ±0.02 -0.99 ± 0.05 41.92 ±0.06 


-2.61 ± 





05 


0.200 < z < 0.300 


-1.20 


42.55 ±0.04 -3.17 ±0.01 -0.88 ± 0.07 42.26 ±0.06 


-2.80 ± 





05 


0.300 < z < 0.377 


-1.20 


42.83 ±0.03 -3.44 ±0.01 -0.66 ± 0.06 42.39 ±0.05 


-2.90 ± 





04 


a The redshift rang 


c 0.010 < z < 0.100 covers an atypical under-dense region (Section [2]) 








Combined SHELS and S08 result. The quoted uncertainties are the formal uncertainties of the fit. 


Note. — For each 


redshift 


range, we list the determined Schechter parameters with fixed a — 




1.20 



(left) and a unconstrained (right) and their uncertainties. We also list the parameters for the Schechter fit 
for a pure star forming sample (top) and a sample where the AGNs and composites are included (bottom). 
The pure star forming sample and the sample including composites and AGNS are based on criterion (iv) 
and (Hi), respectively, in Table [1] The quoted uncertainties are calculated in Section [5.11 
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0.010 < z < 0.100 

0.100 < z < 0.200 

0.200 < z < 0.300 

0.300 < z < 0.377 



58 
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40 
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(ergs E 



42 
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) 



Fig. 10. — Ha luminosity functions for several redshift bins for our pure star forming sample. The Schechter 
functions are derived using the STY-method with a fixed to —1.2. The data-points come from 1/V max where 
the uncertainties are Poisson errors. Table [2] lists the parameters. 



The selection of galaxies by apparent R- 
band magnitude does not yield the same sam- 
ple of galaxies obtained when selecting by Ha 
flux/luminosity. Because we determine the Ha 
luminosity function from an i?-selected spectro- 
scopic survey, we must investigate the poten- 
tial systematic effects of the selection criteria 
(-Rtot < 20.3 and fa a > 10 -155 ergs" 1 cm -2 ) 
on the Ha luminosity function. 

From a given luminosity function (a = —1.20, 
log L* (ergs -1 ) = 42.00, and log 0* (Mpc -3 ) = 
—2.75) we construct a sample of galaxies with a 
flux /hq > 10 -15 ' 5 ergs -1 cm -2 . We choose these 
parameters because they are close to our recovered 
parameters from SHELS. Furthermore, we keep 
the parameters constant over our redshift range 
to test whether our selection criteria introduce an 
artificial evolution to the parameters. We assume 
a uniform galaxy distribution in a comoving vol- 
ume. We assign each simulated galaxy an abso- 
lute magnitude using the distribution of absolute 
magnitude as function of Ln a in Figure [TTJ We 
calculate the apparent magnitude using the allo- 



cated redshift and a k + e-correction based on the 
observed distribution as a function of redshift. We 
apply the survey selection criterion of i?tot < 20.3 
and recover the Schechter parameters using the 
STY-method. Figure [H and Table Ogive the re- 
sults. 

When we keep a fixed in fitting the simula- 
tions, there is an artificial trend of increasing L* 
and decreasing with increasing redshift. This 
trend results from selective removal of the fainter 
Ha galaxies. The removal would otherwise re- 
sult in a decreasing (which we also show for a 
unconstrained). Thus, L* and (f>* should be cor- 
rected for the fact that a is kept fixed at a steeper 
value than would be fit. However, the simulated 
trend in L* and <ft* is far smaller (A logi* = 0.33, 
A log (f>* = —0.20) than we determine from the 
observations (Alogi* = 1.11, Alog^* = -0.73). 
Thus, the evolution in L* in Figure [12] is real. 

When we fit the simulations for all three param- 
eters (unconstrained a in Table [3]), the decrease of 
a (Aa = —1.01) with increasing redshift is close 
to that of the observations (Aa = —0.72; Table[2]). 
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Gallego et al. (1995) 
Tresse & Maddox (1998) 
Sullivan et al. (2000) 
Fujita et al. (2003) 
Hippelein et al. (2003) 
Perez-Gonzalez et al. (2003) 
Ly et al. (2007) 
Morioka et al. (2008) 
Shioya et al. (2008) 
Westra & Jones (2008) 
0.010 < z < 0.100 
0.100 < z < 0.200 
0.200 < z < 0.300 
0.300 < z < 0.377 
SHELS + Shioya et al. (2008) 



Fig. 12. — The three Schechter parameters as a function of redshift and look-bac k tim e for SHELS for pure 
star forming galaxies (red, green, blue, and cyan points), SHELS combined with |S08| (magenta point), and 
surveys at similar redshifts that also use Ha as star formation indicator (black points). 



Moreover, the faint-end slope from our lowest red- 
shift bin (0.010 < z < 0.100; where the faint-end 
of the luminosity function is well-sampled) is con- 
sistent with that of the c ombined luminosity func- 
tion of SHELS and [SOI at z ~ 0.24 (Section 
within the uncertainties. Hence, we have no evi- 
dence for evolution of the faint-end slope over the 
redshift range covered by SHELS. The trend ob- 
served with the faint-end slope unconstrained is 
the result of our selection criteria. 

We also notice an artificial trend in L* with 
increasing redshift for a constant luminosity func- 
tion (although smaller than with a constrained) 
opposite to the trend observed, and opposite to 
the trend we derive fitting for a — —1.20. To com- 
pensate for a shallow faint-end slope and a slight 
decrease in L* with redshift, <jf increases in these 
simulations. 

We do not consider <jf because it only normal- 
izes the luminosity function and does not deter- 
mine the shape of it, unlike a and L* . The nor- 
malization is dependent on the number of galax- 
ies sampled. Because this number is heavily in- 



fluenced by the distribution of galaxies (i.e. the 
large-scale structure, see Figured]), it is not pos- 
sible to say anything meaningful about any trend 
in <fi* even with the area covered by SHELS. 

In summary, there is strong evidence for evolu- 
tion in L* and no evidence for evolution in a over 
0.100 < z < 0.377. 

5.3. The "true" Ha luminosity function 



In Section l5T2l we investigate the influence of our 
selection criteria on an assumed luminosity func- 
tion. We can extend this application to determine 
the "true" Ha luminosity function. 

We construct a sample of galaxies with a flux 
/ho > 10 -15 ' 5 ergs -1 cm~ 2 for a grid of given val- 
ues of a and L* . We constrain <fi* by the number 
of observed galaxies in each redshift bin. These 
choices are our input Schechter parameters. We 
apply our magnitude selection of i?tot < 20.3. 
Then we determine the output: the parameters 
one would recover using the STY-method. We 
take the median of the recovered Schechter pa- 
rameters as our final output Schechter parameters. 
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Table 3 

Recovered Schechter parameters for simulated Ha luminosity functions. 







fixed a 






unconstrained a 




redshift range 


Ot 


logL* 




o. 


logL* 


logr 


input parameters 


— 1.20 


42.00 


-2.75 


— 1.20 


42.00 


-2.75 


0.010 < z < 0.100 


— 1.20 


42.01 ±0.07 


-2.75 ±0.02 


-1.18 ±0.04 


41.98 ±0.10 


-2.71 ±0.08 


0.100 < z < 0.200 


-1.20 


42.06 ±0.03 


-2.77 ±0.01 


-1.07 ±0.03 


41.93 ±0.04 


-2.60 ±0.04 


0.200 < z < 0.300 


-1.20 


42.18 ±0.02 


-2.84 ±0.01 


-0.73 ±0.05 


41.85 ±0.03 


-2.43 ±0.03 


0.300 < z < 0.377 


-1.20 


42.33 ±0.02 


-2.95 ±0.01 


-0.17 ±0.07 


41.81 ±0.03 


-2.38 ±0.01 




Fig. 13. — Histogram of the Ha luminosity function parameters determined from simulations to test the 
influence of the i?-band selection criteria on the Ha luminosity function. Each histogram (red, green, blue, 
cyan) indicates the recovered parameters for a different redshift bin. Vertical dotted lines indicate the input 
parameters of the luminosity function to the simulations (a = —1.20, logL* = 42.00, and log cj)* = —2.75). 
We indicate the results with a fixed at —1.20 (left) and with a as a free parameter (right). Table [3] also 
shows the results. 



With these final output parameters we determine 
the likelihood for our observations (Figure [TU|) . 
We also show the input parameters and the confi- 
dence intervals from the likelihood-determination 
for each redshift bin (Figure [14]) . 

Again, we find a significant evolution of L*. L* 
increased towards higher redshifts, regardless of 
the inclusion of the lowest redshift results. Fur- 
thermore, there is no significant evolution in the 
faint end slope of the intrinsic Ha luminosity func- 
tion. These results confirm the findings in the Sec- 
tion [521 The evolution in L* is real and there is 
no evidence for evolution in a. 

Given the input Ha luminosity function, we 
can calculate the selection function for each red- 



shift bin. We define the selection function as the 
ratio of the measured data points in Figure \TU\ 
and the intrinsic or "true" Ha luminosity func- 
tion. The selection function measures the effect 
of our R < 20.3 selection criterion. We show the 
selection functions in Figure Qjj] 

At z ~ 0.24 we also consider the data of |S08| 
These data should be complete over the luminosity 
range covered by SHELS. We thus assume the data 
from the narrowband survey as the intrinsic Ha 
luminosity function and take the ratio between the 
IS08l and the SHELS data as an estimate of the 
selection function fFigure [T5l magenta long-dashed 
line) . For consistency with the other redshift bins 
of SHELS we remove the OEW Ha+ r Nn j > 12 A 
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Fig. 14. — The luminosity function parameters for the input Ha luminosity function used to determine the 
observed Ha luminosity functions as shown in Figure [TTJl The contours show the 68.3, 95.4, and 99.7% 
confidence intervals based on the fit of the output Ha luminosity function. 



constraint in this calculation (solid diamonds in 
Figure U). 

Th e sel ection function computed at z ~ 0.24 
using |S08| should lie on top of the SHELS selec- 
tion function at 0.200 < z < 0.300, but it does 
not. Thus either SHELS underestimates-or |S08| 
overestimates-the number of faint Ha galaxies. 
Either the i?-band magnitud e vs. Ha-luminosity 
relation between SHELS and |S08| must be signifi- 
cantly different, or there is another selection effect 
not yet considered. Figure [5] rules out a different 
i?-Ha relation. This figure shows that the two sur- 
veys clearly overlap and do not have a significantly 
different i?-Ha relation. 

A selection effect that removes some galaxies 
from the SHELS sample is the OEW HQ+ ^ n j > 

12 A criterion. Figure [9] shows the effect of this cri- 
terion; it slightly increases the number of galaxies 
at the faint-end of the SHELS luminosity function. 
This bias is, however, insuffi cient to explain the 
differences in the SHELS- and lS08l -based selection 
functions. 



The narrowband survey of S08 may overesti- 
mate the number of fainter Ha galaxies. Even 
though consistent within the uncertainties, the 
faint-end slope of the combined luminosity func- 
tion at z ~ 0.24 is somewhat steeper (a ~ —1.4) 
than at our lowest redshift bin (a ~ —1.2) causing 
a very steep selection function at z ~ 0.24. 

We can determine the re dshift of several 
narrowband-survey candidates of lSOSl using zCOS- 
MOS DR2. Figure Hoi shows the redshifts of the 
galaxies from S08J with confirmed redshifts near 
z ~ 0.24. Several galaxies have redshifts outside 
the wavelength range where the NB816 filter is 
sensitive to Ha at z ~ 0.24 (red line). About 25 % 
of the candidates with spectroscopy are at a lower 
redshift z ~ 0.21. This redshift correspond to the 
wavelength range where the NB816 filter is sensi- 
tive to the [Sn] AA6733,6718 doublet (blue lines). 
These galaxies belong to an overdensity in the 
large-scale structure at z ~ 0.22 (dotted histogram 
in Figure [TBI and solid histogram in Figured]). 

Figure [T7] shows the fraction of galaxies with 
redshifts corresponding to [S n] or Ha as a func- 
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Fig. 15. — Logarithm of the selection function for 
each of the redshift bins sampled by SHELS as a 
function of Ha with respect to L* at the respective 
redshift bin. 



tion of the lSOa Ha luminosity. The figure suggests 
that the fraction of [S n] galaxies increases towards 
fainter luminosities. This effect c ould produce an 
excess of faint Ha galaxies in the IS08l survey and 
thus could explain the difference in the selection 
functions implied by the SHELS simulations and 
the comparisons of SHELS and IS08l 

Color-color selections are sufficient to remove 
contaminating galaxies from the narrowband sur- 
vey at higher redshifts. The contaminants include 
H/3 and [O m] at z ~ 0.6 - 0.7, [On] at z ~ 1.2 
and Lya at z ~ 5.7 for a narrowband survey a t 
~ 8150 A (e.g. iFuiita et alj|2003t ILv et al.ll2007l ). 
However, it is impossible to disti nguish Ha galax- 
ies fr om [Sn] galaxies by color ( Westra fc Joned 
2008) . The contamination is survey dependent be- 
cause it depen ds on the details of the large-scale 
structure. The lS08l survey is a case where there is 
a peak in the redshift distribution exactly where 
the narrowband survey is sensitive to [S II] . 



5.4. Volume dependence 

There is a large spread in the parameters de- 
termined from different surveys around z ~ 0.24 
(Figure I12p accompanied by very large uncer- 
tainties. All of these surveys (IFuiita et al. _ 20031 



Hippelein et all2003tlLv et al.ll2007t I Westra fc Jonesl 



2008) use a single or multip le narrowband filters 
950 □' 



S08I uses 5540 □'. 
1 x 10 4 Mpc 3 ; 



over ~ 300 

cal volumes are 0.5 - 1 x 10 4 Mpc 3 ; ISOS 
3 x 10 4 Mpc 3 . The smaller volumes are not large 



Typi- 
covers 



Fig. 16.— Redshifts from zCOSMOS DR2 of Ha 
candidates of IS08l {solid histograms) and all galax- 
ies [dotted histogram), and the relative NB816 fil- 
ter transmission curve for Ha [red li ne) a nd both 
[S n] lines ( blue lines) . Note the |S08| galaxies 
within the [S n] sensitive redshift range. 



enough to constrain the bright end of the lumi- 
nosity function. We discussed IS08l in detail in 
Section gl 

To examine the impact of small volumes, we 
split SHELS into 16 separate pieces to match the 
area (~ 0.25 D°) of typical narrowband surveys 
that probe redshift ~ 0.24. Table 2] gives the me- 
dian recovered parameters and the inter-quartile 
range. 

For a = —1.20 the recovered parameters are 
almost identical to those of the entire field. The 
inter-quartile range is large, even when compared 
to the uncertainties in Table [H If we combine the 
16 "surveys" , we would have to increase the uncer- 
tainties in Table [21 because of the smaller number 
of galaxies, i.e. an increase in shot-noise. This 
uncertainty easily explains the scatter of the pa- 
rameters observed at z ~ 0.24. It underscores 
the need for large- volume surveys to constrain the 
bright end of the luminosity function. 

To constrain a it is more important to have a 
deep survey and to span a large range of luminosi- 
ties r ather than to c over a large area. The data 
from ILv et al. (12007 ) demonstrate this point . As 
discussed in S08 . the data-points from ISOSl and 
Lv et al.l are quite similar at the fainter luminosi- 
ties, both in sl ope and in amplitude. Thus, the 
survey area of ILv et all i.e. ~ 0.25 can be 
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41 .0 

Log Ha -> 



Fig. 17. — Fraction of candidates with a red- 
shift from zCOSMOS DR2 corresponding to Ha 
(z ~ 0.24; red histogram) or [S n] (z ~ 0.21; blue 
histogram) as a function of Ha lumi nosi ty calcu- 
lated from the narrowband survey of IS08l . 

large enough to constrain a, and a only. Their 
area (volume) is too small to determine t he b right 
end of the luminosity function (see also |S08| ) be- 
cause they do not observe enough of the rare most- 
luminous galaxies. 

To estimate the area required to constrain L* 
and (jf , we simulate many observed galaxies given 
a specific luminosity function at 0.233 < z < 0.251 
for different sized areas. We fit the parameters 
with fixed a — —1.20. On average, the param- 
eters are very well recovered. However, for the 
smaller areas the spread in the recovered parame- 
ters is large. We show the 1 a spread around the 
mean, the median, the inter-quartile range, and 
minimum and maximum values of the recovered 
values for each area in Figure [18] 

If we assume that 10 % is an acceptable un- 
certainty for a parameter (~ 0.04 in dex), then 
the survey area r equire d is ~ 3D ° . Su rveys 



like iFuiita et all (l2003r) . iLv et all (|2007t ) and 



Westra fc Jones! (|2008l ) at z ~ 0.24 are thus not 
large enough to con strai n the bright end of the 
luminosity function. |S08| is a factor of two shy of 



this area; SHELS is larger. 

Hence, combining the S08 and SHELS data 
(Section 14. 2|) is an excellent way to constrain the 
faint and bright end of the luminosity function si- 
multaneously. 



Log dQ (degree') -> 

Fig. 18. — Box and whisker plot for the simulated 
surveys as a function of area. The surveys cover 
0.233 < z < 0.251 and have a limiting flux of 
10 -15 5 ergs -1 cm -2 . The gray box indicates the 
1 a around the mean, the dash indicates the me- 
dian, the boxes indicates the inter-quartile range, 
and the whiskers indicate minimum and maximum 
values of recovered Schechter parameters (using 
the STY-method) from the simulations. 

6. Star formation density 

We determine the star formation density (p in 
M Q yr -1 Mpc -3 ) from the integrated Ha luminos- 
ity density for each redshift range. We use the 
conversio n from Ha lumin osity to star formation 
rate from Kennicuttl (1998) for Case B recombina- 
tion and T e = 10 4 K 



SFR = 7.9 x 10 -42 L 



Ha, 



(J) 



where SFR in M Q yr -1 and Lhq in ergs -1 . We 
determine the Ha luminosity density for L > L\\ m 
from the parameters of the Schechter function us- 



ing 



*L*r(a + 2,_^), 



(8) 



where T is the incomplete gamma function. 

The choice Lu m — affects the integrated lumi- 
nosity densitjfl Figure [19] shows the logarithm of 
the total luminosity density evaluated at (log -^p- , 
a) divided by the luminosity density at (log -^f 1 , 
a) = (-1, -1.35). We can thus determine the effect 
of using different limiting luminosities on the total 



6 Aim = reduces Eq (5} to C - 
the complete gamma function. 



t>*L*T(a + 2), where T is 
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Table 4 

Median, upper and lower quartile range of Schechter parameters for 0.25 D° subsets 

using SHELS. 











fixed a 




unconstrained 


a 


redshift 


range 


oc 


logi* 


log>* 


Q 


logL* 


logr 


0.010 < 
0.100 < 
0.200 < 
0.300 < 
0.233 < 


z 
z 
z 
z 
z 


< 0.100 

< 0.200 

< 0.300 

< 0.377 

< 0.251 


— 1 .20 

— 1 .20 
-1.20 
-1.20 
-1.20 


41.60+^if 
42 Ofi +01 ° 

* z - uu -0.09 

42.53t° ^ 
42.82±^ 
42.40_g__ 


-2.79±°;S! 

-2.061SS 

q 90+O.O8 
°- z °-0.12 

-3 5fi+ 05 

°- OD -0.09 

o OO+0.07 

J - zo -0.09 


1 1a +0.l9 

-0 fi2 +01 ° 

u - DZ -0.32 

-o.6ii° ;iI 

-o.48t8:i2 

-0.22t°i 4 


41 c;o+0.23 

41 fi3+ 018 
41.DO_o.22 

42.07t°; 2 2 i 

42.27t° : 15 
41.9ll°oi 7 


O 7R+0.26 
z - ' D -0.34 
9 r 7 +0.28 
z ' ' -0.05 
9 7 r+0.07 
z - ' °-0.07 

-2-791^0 



luminosity density. For example, using L\ nn = 
for a = —1.35 rather than L_ m = 0.1L* gives a 
difference of io 012 ~ 000 = 1.31, i.e. an increase 
of 30%. These effects are obviously more severe 
for steeper values of a. When comparing sur- 
veys of different depths, one needs to be careful 
about extrapolations of the Ha luminosity func- 
tion (Schechter function) to very low star forma- 
tion rates, especially for steep a. 

Table [5] lists the star formation densities and 
uncertainties for SHELS down to the luminosity 
limit of the appropriate redshift bin. We also show 
the star formation densities down to log L\ im = 
40.00 corresponding to a star formation rate of 
0.079 Mq yr _1 for comparison with other surveys 
(Figure |20|). We choose this value for all surveys 
because most surveys either reach this star forma- 
tion rate, or the required extrapolation is modest. 
The solid symbols in Figure represent surveys 
with star formation densities derived from the Ha 
line; the open symbols come from either the [On] 
or [O in] line. Figure [201 shows that the star forma- 
tion density for other surveys at 0.200 < z < 0.300 
is consistent with the star formation density de- 
termined from the combined luminosity function 
of SHELS and [S0l 

Figure [20] also shows a clear increase in the star 
formation density with increasing redshift. How- 
ever, our lowest redshift point (0.010 < z < 0.100) 
lies below surveys at similar redshifts. This un- 
derestimate occurs because our field was selected 
against low redshift clusters. This survey is thus 
an underdense region at low redshifts and the star 
formation density is probably correspondingly un- 
derestimated. 



Because we use the integrated Schechter func- 
tion to determine the star formation density, the 
arguments in Section 15.41 for the Schechter pa- 
rameters L* and <fi*, are valid for the star for- 
mation density. The median, upper and lower 
quartile range for the star formation density for 
0.25 D° subsets are in Table __ Again, the recov- 
ered star formation density of a 0.25 D° subset is 
almost identical to that of the entire field, but the 
standard deviation in the star formation density 
is very large (almost a factor of 2). The large 
uncertainty mainly results from the scatter in L* 
and (jf . Again, combined with the increased un- 
certainties resulting from increased shot-noise, the 
spread in star formation densities at narrow red- 
shift slices can easily be explained by sampling a 
volume that is too small. 

7. Physical properties of star forming 
galaxies 

7.1. Stellar population age 

In star forming galaxies the Ha emission origi- 
nates from gas surrounding the young stars. The 
spectrum from an actively star forming galaxy is 
dominated by the light emitted by these young 
stars. Figure [21] shows Z3„4000, the ratio of the 
continuum red- and bluewards of the H + K break 
and an indicator of the age of the stella r pop- 
ulation ( Bal ogh et al.l [1993 : iBruzual 1983 ). as a 



function of the Ha luminosity for pure star form- 
ing galaxies. A low D„4000 (D„4000 < 1.44; 
D. F. Woods et al. 2010; in preparation) indi- 
cates a young stellar population. The majority of 
the Ha emitting galaxies contain a young stellar 
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Table 5 
Star formation density. 









logp 


log p with log Lii m = 40 


redshift range 


logi lim a 


fixed a 


unconstrained a 


fixed a unconstrained a 



0.010 < z < 0.100 b 37.84 

0.100 < z < 0.200 39.89 

0.200 < z < 0.300 40.55 

0.300 < z < 0.377 40.95 
0.233 < z < 0.251 c 



2.18±0.10 -2.19±0.17 

1.92 ±0.09 -1.92 ±0.12 

1.82 ±0.05 -1.81 ±0.10 

1.81 ±0.03 -1.82 ±0.08 



-2.20 ±0.10 -2.21 ±0.17 
-1.93 ±0.09 -1.92 ±0.12 
-1.81 ±0.05 -1.81 ±0.10 
-1.80 ±0.03 -1.81 ±0.08 
-1.86 ±0.13 



a Li im = AwD 2 L (zi ow ) 

b The redshift range 0.010 < z < 0.100 covers an atypical under-dense region (Section^. 



Combined SHELS and S08 result. 



Note. — We use the Schechter parameters determined for the pure star forming galaxies in Tabled 
We calculate the uncertainties using standard uncertainty propagation for Eq. ([8]) and the uncertainties 
in Table [2 P is in M yr _1 Mpc~ 3 . 



Table 6 

Median, upper and lower quartile range of the star formation density for 0.25 U° subsets 

using SHELS. 



logp 

redshift range fixed a unconstrained a 



0.010 < z < 0.100 
0.100 < z < 0.200 
0.200 < z < 0.300 
0.300 < z < 0.377 
0.233 < z < 0.251 



-2.25 
-1.96 



+0.12 
-0.23 
+0.13 
0.15 

i.84t8;i« 



+0.12 
-0.22 
0.25 



-2.30 
-1-97_ . 15 

1 oc+0.16 

J^-^-o.is 

+0.03 
-0.19 
+0.23 

-0.26 



-1.82 

-1.96 



Note. — The values for the star formation density 
are integrated down to logLi im = 40. 
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Fig. 19. — Logarithm of the total luminosity density evaluated at (log-^™ 1 , a) divided by the luminosity 
density at (log-^f-, a) — (-1, -1.35) (white filled circle). The black lines indicate negative values. 
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Fig. 21. — D n 4000 as a function of Ha luminos- 
ity for each of the four redshift bins for pure star 
forming Ha emitting galaxies. 

population. 

7.2. Galaxy-galaxy interaction 



Sobral et all (<2009h find that the fraction of 



mergers rises with increasing luminosity particu- 
larly around L* . Some of the SHELS galaxies are 



quite luminous in Ha indicating they are under- 
going a starburst. iBarton et al. find that 
a close pass of t wo galaxies ca n initiate a star- 



burst. Following ISobral et all we examine the 



SHELS data to look for evidence of the impact of 
interactions on the Ha luminosity function. Thus, 
we focus on galaxies that may have (or may have 
had) a recent encounter with another galaxy. We 
determine whether each galaxy has an apparently 
nearby "neighbor". 

A galaxy has a neighbor when the velocity 
difference (corrected for redshift) between the 
two galaxies is < 500 km s -1 , and their pro- 
jected separation is < 100 kpc. These values 
are a standard definition of galaxy pairs (e.g 



Barton et all 2000; Patton et al. 2000: Lin et al 



2004; Woods fc Gelled l2007t IPark fc Choil 12009). 



We include the somewhat deeper SHELS catalog 
to look for neighboring galaxies (see Section [2]) . 
This catalog contains spectra of galaxies with 
magnitudes 20.3 < R < 20.6 where the spec- 
troscopy is 52 % complete. The fraction of all our 
pure star forming galaxies that have a neighbor is 
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Fig. 20. — Star formation density as a function of look-back time and redshift for SHELS (red, green, blue, 
and cyan large solid circles) compared with other surveys using the Ha line (solid symbols), or either [ On] 
or [O in] lines (open symbols) as star formation indicator. We also indicate the combined SHELS and IS08l 
point (solid large magenta circle). We calculate the star formation density using the Schechter parameters of 
each survey to a limiting star formation rate of 0.079 M Q yr _1 (corresponds to Lhq = 10 40 ergs^ 1 ) to reduce 
the systematic uncertainty from extrapolation to L = Oergs . 



15.3% (547 out of 35650 

Figure [5U shows the fraction of galaxies with 
one or more neighbors as a function of Ha lumi- 
nosity for the lowest three redshift bins (colored 
histograms) and for the three redshift ranges com- 
bined (thick black histogram). The fraction is al- 
ways a lower limit; deeper spectroscopy might re- 
veal only more neighbors of a galaxy, never fewer. 

It is striking that the fraction of galaxies with 
neighbors increases around L* (and for the lowest 
redshift bin towards the lowest Ha luminosities). 
This result agrees with a rise in the fraction of 
mergers with incre asing Ha luminosity found by 
ISobral et all (2009). The interesting question for 
galaxy evolution is whether the location of the in- 



7 If we decrease our projected separation criterion to 50 kpc, 
the fraction drops with a factor of ~2 to 7.4 % (265 out of 
3565). The fractions in Figure l22l scale with roughly the 
same factor within the uncertainties. Our conclusions are 
not affected by the choice of projected separation. 



crease determines L* , or whether L* determines 
the location of the increase. 

To investigate this behavior further, we investi- 
gate the magnitude difference between the galaxy 
and its neighbor(s). Figure [23] shows the Ha lu- 
minosity as a function of the magnitude difference 
Am R = i?gaiaxy - ^neighboJE We also indicate the 
demarcation b etween minor and maj or pairs, i.e. 
2(e.| 



Woods fc Gell cr 2007). 



\Am\ 

Luminous Ha galaxies with neighbors tend to 
be mostly part of a major pair, and to a lesser 
extent the more luminous galaxy of a minor pair; 
faint Ha galaxies with neighbors can be part of 
a major or minor pair. However, when faint Ha 
galaxies are part of a minor pair, they tend to be 
the fainter (smaller) galaxy. The behavior in Fig- 
ure[23]is consistent with the picture of interaction- 



3 For neighbors in the iJtot > 20.3 catalog we assumed R - 
20.3. Thus, Amu = ^galaxy — 20.3 for these galaxies. 
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Fig. 22. — Fraction of pure star forming galaxies with one or more neighbors for each redshift bin (colored 
histograms) as a function of redshift and the fraction for 0.01 < z < 0.30 (black thick histogram). Colored 
triangles indicate L* (determined with a = —1.20) for each redshift range. 



induced star formation. The increase in the frac- 
tion around L* implies that galaxy-galaxy inter- 
actions are important for the increase of the Ha 
luminosity in these galaxies. 

8. Summary and conclusion 

We use the Smithsonian Hectospec Lensing 
Survey (SHELS) to study Ha emitting galaxies. 
SHELS is complete to R to t = 20.3 over a large 
4D° area. This area yields a large enough vol- 
ume to study the bright end of the Ha luminosity 
function as a function of redshift. 

We determine the Ha flux and attenuation from 
the SHELS spectroscopy. We also identify galaxies 
that host AGNs or are composites. 

We combine the strengths of two surveys, the 
breadth of SHELS (to constrain the bright-end 
of the luminosity fun ction ) and the depth of the 
narrowband survey of IS08I (to determine the faint 
end slope of the luminosity function), to deter- 
mine a well-constrained Ha luminosity function at 
z ~ 0.24. A narrowband survey goes deep over a 



limited field of view to cover the faint end of the 
luminosity function. A broadband selected spec- 
troscopic survey can easily cover a larger volume 
to probe the bright end of the luminosity function. 
The resu lting Schechter parameters are consistent 
with lSO^ within their uncertainties. 

We determine the Ha luminosity function from 
SHELS for four redshift intervals over 0.010 < z < 
0.377. The lowest redshift interval (0.010 < z < 
0.100) covers an atypical underdense region due 
to field selection. The characteristic luminosity L* 
increases as a function of redshift (AlogL* =0.84 
over 0.100 < z < 0.377). 

The star formation density also increases with 
increasing redshift (Alogp = 0.11 over 0.010 < 
z < 0.377). The star formation rate from the 
combined luminosity function of SHELS and S08 
is consistent with that of SHELS alone at 0.200 < 
z < 0.300. 

The fraction of galaxies with neighbors in- 
creases by a factor of 2 — 5 around L* for the 
most luminous star forming galaxi es at each red- 
shift, similar to Sobral et a l. (2009). The fraction 
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Fig. 23. — Magnitude difference between the galaxy and its neighbor as a function of Ha luminosity for pure 
star forming Ha emitting galaxies. The solid lines show |Arrift| = 2, the demarcation between major and 
minor interactions. Galaxies above the dotted line are fainter than their neighbor, galaxies below are more 
luminous. 



appears to also increase towards fainter Ha lumi- 
nosity as a result of interactions in minor pairs. 
We conclude that triggered star formation is im- 
portant for both the highest and lowest luminosity 
Ha galaxies. 

The future of surveys for star forming galaxies 
is a combination of a large-area spectroscopic sur- 
vey combined with very deep narrowband imag- 
ing. However, the narrowband imaging requires 
extensive test spectroscopy because the impact of 
large-scale structure with respect to the filter re- 
sponse is unknown a priori. The combination of 
methods can constrain and remove the scatter in 
the star formation density as a function of redshift. 
The combination also allows a secure determina- 
tion of the shape of the luminosity function over 
a large luminosity range. 
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